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We study the interplay between superconductivity and antiferromagnetic spin fluctuations in 
electron-doped cobaltate Na^OoC^ • J/H2O based on the kinetic energy driven superconductivity 
mechanism. We show that the superconducting state is governed by both charge carrier pairing 
and quasiparticle coherent, and displays a common dome-shaped phase diagram in agreement with 
experimental results. By calculating the dynamical spin structure factor, we theoretically find that 
the magnetic excitation shows a commensurate resonance peak, which locates at antiferromagnetic 
wave vector Q(^, for a broad range of low energies, then evolves outward into six incommen- 
surate magnetic scattering peaks with increasing energy. Such commensurate-incommensurate spin 
resonance excitation should be measured by the inelastic neutron scattering technique (INS). Our 
present results strongly suggest that magnetic resonance can indeed be one of the fundamental 
features in doped Mott insulators. 
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I. INTRODUCTION 

In the past decade, superconductivity in sodium cobal- 
tate Na^CoC^ • yH 2 (NCO) has been intensively stud- 
ied in that it shares many similarities with high-T c 
cuprateiPHll The doped cobaltate has a lamellar struc- 
ture consisting of two-dimensional (2D) triangular C0O2 
layers separated by a thick insulating layer of Na + ions 
and H2O molecules^. Although the T c is a little bit 
lower, it displays the same kind of doping-controlled 
phase diagram that is observed in the doped cuprates^l 
In particular, the electron-doped NCO may be another 
unconventional superconductor coming from the fam- 
ily of doped Mott insulators. Although the ferromag- 
netic correlations is present in Na x Co02 for large elec- 
tron doping regimes^, the antiferromagnetic (AF) short- 
range spin correlations in Na:i;Co02 and NCO at low 
doping levels has been observed from quadrupolar res- 
onance and t herm opower as well as other experimental 
measurements^^!!!] Therefore, superconductivity medi- 
ated by the AF spin fluctuations in the doped cobaltate 
suggests the importance of the strong electron correla- 
tion as in the doped cuprates^l Understanding its su- 
perconducting (SC) mechanism could provide important 
insight into the microscopic origin of the unconventional 
superconductivity. 

Recently, the kinetic energy driven SC mechanism ap- 
plied into the doped cuprates base d on the charge-spin 
separation fermion-spin theorjl^USI nas been successfully 
studied, where the dressed charge carriers interact oc- 
curring directly through the kinetic energy by exchanging 
spin excitations, leading to a net attractive force between 
dressed charge carriers, then the electron Cooper pairs 
originating from the dressed charge carriers pairing state 
are due to the charge-spin recombination, their conden- 
sation reveals the SC ground state, and the SC transi- 
tion temperature is controlled by both charge carrier gap 



function and quasiparticle coherent weight. Since there 
is a remarkable resemblance in the normal and SC -sta te 
properties between NCO and the doped cupratePH^ as 
mentioned above, and since the strong electron correla- 
tion is common for both materials, it thus could be ex- 
pected that two systems may have similar underlying SC 
mechanism, i.e., superconductivity in the electron-doped 
cobaltate is also driven by the kinetic energy. Taking it 
into consideration, we in this paper study the unconven- 
tional superconductivity of sodium cobaltate NCO, and 
find that the SC state is governed by both charge carrier 
pairing and quasiparticle coherent, and the resulting SC 
transition temperature is suppressed to a lower temper- 
ature due to the strong magnetic frustration, which is 
in agreement with experimental measurements^. Then 
we theoretically analyze the feedback effect of super- 
conductivity with d\ + «c?2-wave pairing symmetry on 
the magnetic excitations, and find a commensurate res- 
onance peak locating at antiferromagnetic wave vector 
Q(^, ^7=) for a broad range of low energies, then evolv- 
ing outward into six incommensurate magnetic scatter- 
ing peaks with increasing energy, which is very similar 
to the case of electron-doped cuprates^^, and suggests 
that magnetic resonance could be intimately related to 
superconductivity. We propose that the inelastic neutron 
scattering experiment can verify the magnetic resonance 
excitations in electron-doped cobaltates NCO after the 
overcome of the difficulty in growing its high-quality SC 
sample. 



II. THE t-J MODEL AND KINETIC ENERGY 
DRIVEN SC THEORY 

In electron-doped cobaltate, the characteristic feature 
is the presence of the 2D C0O2 plane^, of which the essen- 
tial physics is contained in the t-J model on a triangular 
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lattictPMI as 



where C\ a (Ci a ) is the fermion creation (annihilation) op- 
erator subject to the no-double-occupancy constraint, 
Si = CjaCi/2 is the spin operator with a = (cr x , a y ,a z ) 
as the Pauli matrices, /i is the chemical potential. Note 
here that we work in the hole representation, so that a 
hole in the above model represents a physical electron- 
double-occupancy in electron-doped cobaltate. 

We apply the charge-spin separation fermion-spin 
theorySS within which the hole operators can be de- 
coupled as C*f = a\^S^ and = a\^S^ , where the 
spinful fermion operator at a — e _l * <CT a.j describes the 
charge degree of freedom together with some effects of 
the spin configuration rearrangements, while Si describes 
the spin degree of freedom, then the single occupancy lo- 
cal constraint J2a C L C iT = S t a it a lt S i~ +S ^ a il a ll S i^ = 

aia\(S^S~ + S~S^) = 1 — ajaj < 1 is satisfied in ana- 
lytical calculations. Then the low-e nergy behavior of the 
t-J model can be expressed a d 15 * 16 !, 



t') = ((al t (t);al(t'))) aS 

g(k) = gM (A:) + g <°> (fc) [£<*> (k)g(k) - E< a) (-fc)Qft (fc)] , ( 4 ) 
9ft (fc) = g^(-k)[4 a \-k)^(-k) + E; (a) (-%(fc)],(5) 

respectively, where the four vector notation k = (icJ n , k), 
the mean-held (MF) dressed fermion Green's func- 
tion g(°' l (k) — iuj n — £k, with the MF dressed 
charge carrier excitation spectrum £k = Ztx^k — /A 
Z is the number of the nearest neighbor sites, 7k = 
(1/Z)^e 4k ^ = [cosfc, + 2cos(A: :c /2)cos(v / 3fc y /2)]/3, 
while the self-energies have been derived as, 
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with J c ff = (1— x) J, and the doping level x = (a^aicr) = 
(a\ai). In this case, the magnetic energy (J) term in t- 
J model is only to form an adequate dressed spin con- 
figuration, while the kinetic energy (t) term has been 
transferred as the dressed charge carrier-spin interac- 
tion, which dominates the essential physics^. This 
dressed charge carrier-spin interaction is quite strong, 
and induces the dressed charge carrier pairing state (then 
the electron pairing state and superconductivity) by ex- 
changing spin excitations in a higher power of the elec- 
tron doping level as in the doped cupares^. Then the or- 
der parameter for the electron Cooper pair without AF 
long range order (AFLRO) can be expressed as, 



(c! t cl-ciq t ) = -{s?sr)A a 
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with the dressed charge carrier pairing order parameter 
A a = (a^a,^ — aj-fai±). This shows that the dressed 
charge carrier pairs move freely in the background of 
the disordered spin liquid state, and then the phys- 
ical properties of SC state are essentially determined 
by the dressed charge carrier pairing state. Based on 
the Eliashberg's strong coupling theory^, we obtain 
the self-consistent equations in terms of the method of 
equation of motion that is satisfied by the full dressed 
charge carrier diagonal and off-diagonal Green's func- 
tions g(i-j,t-t') = ((a ia {t);a] a {t'))) and &(i-j,t- 



where the MF dressed spin Green's function D^(i — 
j,t — t') = ((5j + (t); S~(t')))o- Since the pairing force and 



dressed fermion gap function have been incorporated into 
the self-energy E^^fc), then it is an effective dressed 
charge carrier gap function. While self-energy T,^\k) 
renormalizes the MF dressed charge carrier spectrum, 
so that it describes quasiparticle coherence and domi- 
nates the charge transport of the materials^. Moreover, 
T,^\k) is an even function of iui ni while the self-energy 
E^ (fc) is not, which can be split into its symmetric and 



antisymmetric parts as E^ (A;) = Y,f e '(k) + iuj n Y^J{k) 
with E^" (k) and E^ (k) being even functions of iui n . 
Now we define the charge carrier quasiparticle coherent 
weight Z F \k) = l-i;[ a J(k). 

We only study the static limit of the effective dressed 
charge carrier gap function and quasiparticle coherent 

weight, i.e., E^(fc) = A Q (k), and Z F x (k) = 1 - E^(k). 
Then the dressed charge carrier diagonal and off-diagonal 
Green's functions are obtained as, 
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with E^k) = E_ k , E 2 (k) = -E kl & = Z F (k)&, 
A* z (k) = Zp(k)A*(k) and the dressed fermion quasi- 
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particle spectrum E k = ^£k"+ | A a z(k) | 2 . Although 
Z F (k) is a function of k, the wave vector dependence 
is unimportant in that everything happens at the fermi 
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surface (FS). In this case we approximate Zp(k) by a 
constant Zp near the F£P. So far the pairing symmetry 
in NCO is far from reaching a consensus, many experi- 
mental data point towards to spin-singlet Cooper pairing 
and suggest non-s-wave superconductivity without a full 



gap 1 2 * 3 * 7 * 8 -! In particular, it has been argued according 
to the irreducible representations of the triangular lat- 
tice that there are three possible basis functions of even 
parity^, i.e., one s-like function Sk = cosfe^ + cos[(A; x — 
v3fcj/)/2] + cos [(fc^ + y/3ky)/2], and two d-like functions, 
dik = 2cosfcic — cos [(k x - V3k y )/2] - cos[{k x + V3k y )/2] 
and <i 2 k = V3cos{(k x + \fSk y )/2] — v / 3cos[(fc 2 , — y/3k y )/2]. 
However, with the different linear combinations of these 
basis functions, it has been found^l in terms of the 
Gutzwiller approximation scheme and variational Monte 
Carlo simulation that the lowest energy state is the d- 
wave (di -\-1d2) state. Therefore we in the following calcu- 
lations consider both extended s-wave A a s (k) = 

^1 ^ 5k 

and the d-wave case A^(k) = A„ (dik+^k), and then 
the effective dressed fermion gap parameter and quasi- 
particle coherent weight satisfy the following equations, 
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respectively, where r = s,d, F^ (k, q, p) = 
(uj p - w q )[n fl (w q ) - ra B (<j p )][l - 2n F (F k )] + 

F k [n B (w p )n B (-w q ) + n B (w q )n B (-w p )], Ff ; (k,q,p) = 
"K + cj q )[n B (w q ) - ns(-w p )][l - 2n F (F k )] + 

F k [n B (w p )n B (iLj q ) + n B (-CJp)n B (-u;q)], F 2 (1) (q,p) = 
n F (F p _ q+ko )[n B (w q ) - ji b (w p )] - n B (w p )n B (-w q ), 

("2") 

^2 (q.P) = n F (E p _q +ko )[n B (ujp) - n B (ujq)} - 
n B (ujq)n B (-ujp), F$ ; (q,p) = n F (F p _ q+ko )[n B (cj q ) - 
n B (-w p )] + n B {LOp)n B {ujq), and F 2 (4) (q,p) = 
"F(F p _ q+ko )[n B (-o; q ) - n B (ujp)} + n B (-<j p )n B (-Wq)- 
These two equations determine the SC order directly, 
and must be solved simultaneously with other self- 
consistent equations, then all order parameters are 
determined by the self-consistent calculation^. In terms 
of the off-diagonal Green's function (9), we obtain the 
dressed charge carrier pair order parameter as, 



A i T) = ^El7 k T) | 2Z ^tanh[> u -. ,1.2, 



The dressed charge carrier pairing state originating from 
the kinetic energy term by exchanging dressed spin exci- 
tations will also lead to form the electron Cooper pair- 
ing state as mentioned in Eq. (3). For a discussion 
of the physical properties in SC state, we now need 
to calculate the electron off-diagonal Green's function 
rt(i -j,t- t') = «Cj.(t); Cjj. (*')))• In thc framework 
of the charge-spin separation fermion-spin theory^, it is 
a convolution of the dressed spin Green's function (p) 
and off-diagonal dressed charge carrier Green's function 
^s{k), and can be expressed as, 

p ip m 

1 ^24 T i(p-k) Bp 



N 
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P -k 



2w r _ 



(uj p + £p_ k )[ng(M p ) + n F (-F p _ k )] 
(iu n ) 2 - (oj p + F p _ k ) 2 
(up - Ep- k )[n B (Ljp) + n F (F p _ k )] 
(iuj n ) 2 - (w p - F p _ k ) 2 

then the SC gap function is obtained as, 



(13) 



A (T) (k) = -i^rt(k,io;„) 



lE ^( P -k) 1 

2F p _ k 2 1 p kJ 



B 1 

X 2£ C ° th[ 2^ p] ' 



(14) 



F k 



which can be further written as A^ T )(k) = A^ r ^7^. With 
the help of Eqs. (14) and (12), the SC gap parameter can 

be obtained as A^ r ^ = — %Ai T \ then the effective SC gap 
parameter should be A( T ) ~ — ^A^'. As in the doped 
cupratea^, our present theory clearly indicates that there 
is a coexistence of the electron Cooper pair and short- 
range AF correlation in the doped cobaltates, and there- 
fore the short-range AF fluctuation can persist into su- 
perconductivity. The reason is that the AF fluctuation is 
dominated by the scattering of dressed spins^, which has 
been incorporated into the electron off-diagonal Green's 
function (and hence the electron Cooper pair) in terms 
of the dressed spin Green's function. This result is con- 
sistent with some experimental results^^HID 

We perform a numerical calculation for the gap pa- 
rameters, and the results of the dressed charge carrier 
(a) and SC (b) gap parameters in the s-/d-wave symme- 
try as a function of the electron doping x at T = 0.001J 
and t/J = —3.0 are shown in Fig. 1, where both values 
of the dressed charge carrier and SC gap parameters in- 
crease with increasing doping in the underdoped regime, 
and reach the maximal value for a particular doping con- 
centration, then decrease in the overdoped regime. Al- 
though there is a coexistence of the electron Cooper pair 
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FIG. 1: (Color online) The effective dressed charge carrier 
pairing (a) and effective superconducting (b) gap parameters 
in the s-wave symmetry (green line) and d-wave symmetry 
(red line) as a function of doping concentration in T = 0.002J 
and t/J = —3.0. The superconducting transition temperature 
(c) as a function of doping concentration in the s-wave sym- 
metry (green line) and d-wave symmetry (red line). Inset: 
the experimental result taken from Ref. [2]. 



and short-range AF correlation, the value of the SC gap 
parameter is still suppressed by this AF fluctuation. 

The SC transition temperature T c occurs when A*- 7 ") = 
0, which is identical to the dressed charge carrier pair 
transition temperature occurring in the case of A 1 //' 1 = 0. 
The SC transition temperature T c as a function of dop- 
ing x in the s-/d-wave symmetry for t/J — —3.0 is plot- 
ted in Fig. 1 (c) in comparison with the experimental 
data^ taken from NCO (inset). For the s-wave symme- 

(s) 

try, the maximal SC transition temperature Tc occurs 
around a particular doping concentration x f=a 0.21, and 
then decreases in both lower doped and higher doped 
regimes. However, for the d-wave symmetry, the max- 
imal SC transition temperature T^ occurs around the 



optimal doping concentration x opt ~ 0.29, and then de- 
creases in both underdoped and overdoped regimes. It is 
clearly shown that the SC state has the d-wave symme- 
try in a wide range of doping, in qualitative agreement 
with the experimental data 3 supporting the domina- 
tion of d\ + ic?2 pairing symmetry in NCO. Using an 
reasonable estimation value of J ~ 15 mev to 20 mev 
in NCO, the SC transition temp erature in the optimal 
doping is T°P timal « 0.02J « 4K! 14 ' 15 ' 20 ! Since NCO is 
the electron doped Mott insulator on a triangular lattice, 
the system has strong geometrical spin frustration. This 
magnetic frustration also induces the strong charged car- 
rier's quantum fluctuation. In comparison with the case 
in the doped cuprates, we find that the SC transition 
temperature in the doped cobaltates is suppressed heav- 
ily to a lower temperature due to both strong magnetic 
frustration and dressed charge carrier's quantum fluctu- 
ation. 



III. ENERGY DEPENDENT COMMENSURATE 
AND INCOMMENSURATE MAGNETIC 
RESONANCE 



We now turn to discuss the feedback effect of super- 
conductivity with the e?i + ic?2-wave pairing symmetry on 
the magnetic excitations in the doped cobaltates. Among 
various phenomena of the spin excitations in SC state, 
the neutron spin resonance or the magnetic resonance ob- 
served in SC cuprates is undoubtly the most important 
feature, which has been regarded as an essential feature 
in doped Mott insulators. Unfortunately, due to the dif- 
ficulty in growing the high-quality NCO SC sample there 
is no inelastic neutron scattering experiments on this sys- 
tem, so that it will be timely and necessary to study this 
important issue from the point of theoretical view. 

Within the kinetic energy driven SC mechanism, the 
AF fluctuation is dominated by the scattering of dressed 
spins with the full dressed spin Green's function 



D(k,w) = 



L»(°)- 1 (k,w)-S( s )(k,a;) : 



(15) 



where the second order spin self-energy E^ (k, oS) is ob- 
tained from the dressed charge carrier bubble in terms of 
the collective mode in the dressed charge carrier particle- 
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particle channel as, 

^\k,uj) = {Ztf^-^h 

q 

Bq+k4Ai1(p)Ai1(p + q) 



2 2 
jy2 'q+p+k + 7p-kJ 

p,q 



E p E p+q 



^ (1) (k,p,q) 



+ 



(E p - E p+q 
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p+q ~Ep+ ^q+k) 2 

p(3) 



E s W (k,p,q) 



W 2 - (-Ep + E p+q + Wq+k) 
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p+q 
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q+kj 
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where F« (k, p, q) 

Wq+k){lB (W q+ k) [n F (Ep) 

«-F(E p+q )n F (-E p )}, 

1- w q+ k){"-B(wq + k)[riF(Ep +q 



E, 



p+q 



Ei 2) (k,p,q 



E, 



(E P - 
n F (E p+q )] 

= (Ep +q 

- n F (E p )} 

= (£ P 
- n F (E p+q )] 

n F (-E p+q )n F (-E p )}, E s (4) (k,p,q) = (E p 
w q+ k){"-B(w q+k )[n F (-Ep) - n- F (E p+q )] 



n F (E p )n F (-E p+q )}, 



Ei 3) (k,p,q 



F 



p+q 



J q+k ) {n B (w q+k ) [n F (-E- 



E, 



p+q 



n F (Ep_|_q)n F (E p )}. Then the dynamical spin structure 
factor which determines the spin excitations in the SC 
state can be written as, 



S(k,u) =-2[l 



riB(w)]im£>(k, uj) 
2[1 + n B (cj)]E 2 ImI]( s )(k,w) 



[w 2 - cj 2 - E k ReE( s )(k, w)] 2 + [E k Im£( s )(k,u; 

where ImE^k,^) and Re£( s >(k,w) are imaginary and 
real parts of the second order self-energy respectively. It 
is shown that the above dynamical spin structure fac- 
tor has a well-defined resonance character, where S(k, uj) 
exhibits resonance peaks when the incoming neutron en- 
ergy uj is equal to the renormalized spin excitation, i.e., 
uj 2 — u; 2 . — Ek c Re£( s - ) (k c , uj) = for certain critical wave 
vectors k c , then the weight of these peaks is dominated 
by the inverse of the imaginary part of the dressed spin 
self-energy 1/Im£( s )(k c , uj). 

We plot in Fig. 2 the dynamical spin structure fac- 
tor S^k, uj) of the electron-doped cobaltates NCO in the 
(k x , k y ) (here we use the unit of [7r/3, tt/\/3}) plane at the 
optimal doping x opt = 0.29 with temperature T = 0.002 J 
for parameter t/J = —3.0 at energy (a) uj = 0.08 J, and 
(b) uj — 0.20J. We find that a spin resonance peak lo- 
cates at antiferromagnetic wave vector Q(^, at low 
energy, then evolves outward into six incommensurate 
magnetic scattering peaks at high energy. Along one di- 
rection k y (k x — y), we calculate the energy dependence 
of S(k,uj) and find in Fig. 3 that a commensurate spin 
resonance peak locating at antiferromagnetic wave vector 




FIG. 2: (color online) The dynamical spin structure factor 
S(k, cu) in the (k x , k y ) plane as a unit of [ir/3, tt/VS] at x opt = 
0.29 with T = 0.002J for t/J = -3.0 at (a) uj = 0.08J, and 
(b) uj = 0.20J. 



' v§) P ers ists in a broad range of low energies up to 
about 0.16 J «3.2 mev, then splits into two incommensu- 
rate magnetic scattering peaks with increasing energies 
(The energy dependence of S'(k, uj) along the other two 
directions display the same transition feature and don't 
show here) . The present result is very similar to the spin 
resonance excitation of electron-doped cupateJ2H^, and 
differs from the case of hole doped cuprate in which an- 
other incommensurate-commensurate spin reson ance ex- 
citations transition happens at low energies^^H. There- 
fore the commensurate magnetic resonance could be inti- 
mately related to superconductivity and is a common and 
essential feature in doped Mott insulators, other details 
such as the incommensurability and hour-glass dispersion 
found in different cuprates superconductors may not be 
fundamental to superconductivity. 
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the mechanism of kinetic energy driven superconductiv- 
ity. We find that the SC state is governed by both 
charge carrier pairing and quasiparticle coherent, and 
the resulting SC transition temperature is suppressed to 
a lower temperature due to the strong magnetic frus- 
tration, and shows a common doping-dependent dome- 
shaped feature, in agreement with experimental measure- 
ments. Then we theoretically analyze the feedback effect 
of superconductivity with d\ + ic?2-wave pairing symme- 
try on the magnetic excitations, and find a commensu- 
rate resonance peak locating at antiferromagnetic wave 
vector Q(^f, ^7=) for a broad range of low energies, then 
evolving outward into six incommensurate magnetic scat- 
tering peaks with increasing energy, strongly suggesting 
that magnetic resonance could be one of the fundamen- 
tal features in doped Mott insulators. We propose that 
the inelastic neutron scattering experiment can verify the 
magnetic resonance excitations in electron-doped cobal- 
tates NCO after the overcome of the difficulty in growing 
its high-quality SC sample. 



FIG. 3: (color online) The energy dependence of dynamical 
spin structure factor S(k, w) along the direction k y (k x — 4^) 
in first BZ at x opt = 0.29. 

IV. SUMMARY AND DISCUSSIONS 

In summary, we investigate the unconventional super- 
conductivity in electron-doped cobaltates NCO based on 
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